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Abstract 

We study a holographic model with vector condensate by coupling the anti- 
de Sitter gravity to an Abelian gauge field and a charged vector field in (3 + 1) 
dimensional spacetime. In this model there exists a non-minimal coupling of 
the vector field to the gauge field. We find that there is a critical temperature 
below which the charged vector condenses via a second order phase transition. 
The DC conductivity becomes infinite and the AC conductivity develops a gap 
in the condensed phase. We study the effect of a background magnetic field on 
the system. It is found that the background magnetic field can induce the con- 
densate of the vector field even in the case without chemical potential/charge 
density. In the case with non-vanishing charge density, the transition tempera- 
ture raises with the applied magnetic field, and the condensate of the charged 
vector operator forms a vortex lattice structure in the spatial directions per- 
pendicular to the magnetic field. 
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1 Introduction 


The gauge/gravity duality [1, 2, 3] turns out to be a useful and complimentary framework 
to study a strongly coupled system through an appropriate gravity theory living in a 
higher dimensional spacetime. There are two complementary approaches. In the top- 
down approach, one can obtain some low energy effective theories in the bulk as consistent 
truncations of string/M theory by dimensional reduction on some compactified manifolds. 
The advantage of this approach is that one knows the origin of the effective gravity theory 
in the bulk and the details of the dual field theory. One of the well-known examples in this 
category is the duality between the IIB superstring theory on the AdS; x S° and the N = 4 
supersymmetric Yang-Mills theory on the AdS; boundary. On the other hand, in the so- 
called bottom-up approach, the gravity theory in the hulk is usually constructed by some 
physical considerations, according to the AdS/CFT dictionary. The models constructed 
in this way are simple and universal, just like the Landau-Ginzburg theory describing 
superconductivity. The disadvantage is that the details of the dual field theory are not 
very clear. 

A well known example in the second category is the holographic superconductor model [4, 
5]. More specifically, to build a holographic superconductor model, one should first intro- 
duce a U(1) gauge field in the bulk, which corresponds to a global U(1) symmetry in the 
boundary side. To break this U(1) symmetry spontaneously, one needs a charged operator 
condensing at low temperature. Therefore, one includes a charged scalar field in the bulk 
which is dual to the boundary scalar operator. For simplicity, the gauge field and the 
charged scalar can be minimally coupled. This holographic toy model admits black holes 
with scalar hair at low temperatures (superconducting phase), but without scalar hair at 
high temperatures (normal phase). In this way, one mimics the superconductor /conductor 


phase transition in the field theory side. Holographic superconductor models constructed 
in the top-down approach can also be found, for example, in refs. [6, 7, 8, 9]. 

It is well known that superconductivity involves the formation of a quantum condensate 
state by pairing conduction electrons. The pair of electrons, called the Cooper pair, can 
be in a state of either spin singlet with total spin s = 0 or spin triplet with s = 1. The 
order parameter in a superconductor is expressed in terms of the gap function. Due to the 
anti-commuting properties of the electron wave function, the antisymmetric spin-singlet 
state is associated with a symmetric orbital wave function with orbital angular momentum 
l = 0 (s-wave), 2 (d-wave), etc, while the symmetric spin-triplet state is accompanied by 
an antisymmetric orbital wave function with orbital angular momentum l = 1 (p-wave), 
3 (f-wave), etc. Since the condensed field in the holographic setup [4, 5] is a scalar field 
dual to a scalar operator in the field theory side, it is therefore a holographic s-wave 
model. In the holographic setup, a d-wave order parameter is dual to a charged spin 
two field propagating in the bulk [10, 11]. To mimic p-wave superconductor, one can 
consider a vector order parameter which is dual to a vector field in the gravity side. A 
holographic p-wave model [12] was constructed by adding a SU(2) Yang-Mills field into 
the bulk. A U(1) subgroup of the Yang-Mills field is regarded as the gauge group of 
electromagnetism. A gauge boson generated by another SU(2) generator charged under 
this U(1) by the nonlinear coupling of the non-Abelian field is dual to the vector order 
parameter. An alternative holographic realization of p-wave superconductivity emerges 
from the condensation of a 2-form field in the bulk [13]. A holographic chiral ps + ip, 
superconductor was discussed by generalized the SU(2) model introducing a Maxwell field 
and a Chern-Simons term [14]. Note that in principle we can also build a holographic 
p-wave model by introducing a complex vector field charged under a U(1) gauge field in 
the bulk with possible couplings. 

On the other hand, motivated by the possibility to create a very strong magnetic 
field, for example, at RHIC and LHC, the interest to investigate the properties of QCD 
matter in a strong magnetic field has been growing recently. Some new phenomena have 
been revealed such as the chiral magnetic effect, a split between the deconfinement phase 
transition and chiral symmetry restoration phase trasnition, see ref. [15] for a review. An 
interesting new phenomenon is the possibility that the QCD vacuum undergoes a phase 
transition to a new phase with charged p-meson condensed in a sufficiently strong magnetic 
field [16, 17]. This exotic phase is a kind of anisotropic superconducting phase [18]. The 
author in ref. [16] adopted an effective quantum electrodynamics action (DSGS model |19]) 
to discuss the condensate of p-meson. A similar study based on Nambu-Jona-Lasinio 
model [20] can be found in ref. [17]. The uniform magnetic field background is encoded 
in an Abelian gauge field. As a vector boson, the possible condensate of the p-meson in 
a uniform magnetic field can be realized in the holographic framework by introducing a 
charged vector field in the gravity side in the presence of a magnetic field. Based on the 
Sakai-Sugimoto model, the holographic p-meson was studied in the confinement phase 
at zero temperature [21] where it is shown that the effective mass of the p-meson at 
strong magnetic field becomes tachyonic. The influence of the magnetic field on the chiral 
transition temperature and deconfinement transition temperature was discussed in ref. [22]. 
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In the SU(2) model [12], a similar instability triggered by a non-Abelian magnetic field 
has been found in refs. [23, 24, 25], which is reminiscent of the observation that non- 
Abelian magnetic field induces the W-boson condensate exhibiting vortex lattices in flat 
spacetime [26, 27, 28]. 

In this paper we will construct a toy model by introducing a complex vector field 
charged under an Abelian gauge field in the bulk, which is dual to a strongly coupled 
system involving a charged vector operator with a global U(1) symmetry. In this model 
there exists a non-minimal coupling between the vector field and the Abelian gauge field. 
In this bottom-up approach, we do not know the details of the dual field theory clearly. 
However, our setup meets the minimal requirement to construct a holographic p-wave 
superconductor model. This model has also potential to discuss the p-meson condensate 
in the presence of magnetic field. 

We first study the model in the case without a background magnetic field, i.e., turn 
on the temporal component of the gauge field only, which is the usual way of analyzing 
holographic superconductor at a first step [4, 10, 11, 12]. We find a critical temperature 
below which a vector operator acquires a vacuum expectation value breaking the Abelian 
symmetry spontaneously. Furthermore, the condensate of this vector operator picks out 
a special spatial direction, thus the rotational symmetry is also broken in the condensed 
phase. The system undergoes a second order phase transition with the critical exponent 
one half which coincides with the result from mean field theory. This condensed phase 
presents characteristics known from superconductivity, such as an infinite DC conductivity 
and a gap in the optical conductivity. 

An instability induced by a background non-Abelian magnetic field has been reported, 
where the non-Abelian current operators obtain vacuum expectation values resulting in a 
vortex lattice structure [24, 25]. We are interested in how the applied Abelian magnetic 
field influences the instability of our model with Abelian gauge symmetry. For this, we 
turn on a uniform magnetic field B in the bulk, which immerses the condensed phase into 
an external magnetic field. 

We find that the increase of the magnetic field induces the instability of the black hole 
background, which gives rise to a family of condensate induced by the applied magnetic 
field in the dual strongly coupled system. This magnetic field induced instability can 
happen even for the case with vanishing chemical potential or charge density. For finite 
magnetic field, there is a tower of “droplet” solutions [29] in the sense that they are localized 
in a finite region. Further, the condensate shrinks in size as one increases the magnetic 
field, which shares similarity with the Meissner effect. But it is not the exact case with 
conventional superconductivity, where the superconductivity is suppressed by magnetic 
field in the real materials. Of course, it is interesting to notice the fact that it has been 
reported recently that an applied magnetic field may induce superconductivity [30, 31]. In 
addition, the emergence of charged vector operator condensate triggered by the applied 
magnetic field is consistent with the appearance of electromagnetically superconducting 
phase in a strong magnetic field studied from both field theory method [16, 17] and the 
holographic setup [21, 22]. We also manage to construct the vortex lattice solution in 
the condensed phase near the phase transition, which is very reminiscent of the Abrikosov 
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lattice in common type-II superconductors. 

The plan of this paper is as follows. In section 2, we introduce the holographic model 
and deduce the general equations of motion of the system. In section 3, we turn on a 
uniform magnetic field and give details about how to recover the Landau levels in this 
holographic model. Section 4 is devoted to discussing the phase diagram with/without 
chemical potential. The vortex lattice solutions are constructed in Section 5. The con- 
clusion and further discussions are included in section 6. We study, in Appendix A, the 
condensate of the vector field and the phase transition in the case without magnetic field 
by calculating the conductivity to ensure the condensed phase to be a superconducting 
state. 


2 The holographic model 


Let us introduce a charged vector field into the (3+1) dimensional Einstein-Maxwell theory 
with a negative cosmological constant. The full action reads 


1 YV 1 V 
= S= 55 | dev IR +E- FoF" — plap” — mpo" + igy F], (1) 
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where L is the AdS radius that we will set to be unity, k? = 87G is related to the 
gravitational constant in the bulk and m is the mass of the charged vector field p,. The 
strength of U(1) field A, is Fw = VuAv — VvA,. The tensor p,, in (1) is defined by 
Pu = Dupv — Dy py. The covariant derivative D, = V,,—1iqA, with q the charge of vector 
Py. The last interacting term describes the non-minimal coupling of the charged vector 
field p, to the U(1) gauge field A,,. The parameter y characterizes the magnetic moment of 
the vector field p, and is assumed to be non-negative. The form of the action is reminiscent 
of the DSGS model which describes the quantum electrodynamics of the p-meson proposed 
by Djukanovic, Schindler, Gegelia, and Scherer in ref. [19]. Compared to the action of the 
DSGS model, it is easy to find that the action in (1) is just a simple generalization of the 
DSGS model to the anti-de Sitter space. Note that the part of neutral p-meson in ref. [19] 
is neglected here since it is not relevant to our goal in this paper. 

Varying the action (1) with respect to A,, yields the equation of motion for gauge field 


V” Fin = ial oh, — P Puu) + igyV" loet, — phon), (2) 


while a variation of the action (1) with respect to p”t gives the equation of motion for the 
charged vector field 

D” py — Mpu + ig 0" Fin = 0. (3) 
If one takes the limit q — oo keeping gp, and qA, fixed, the back reaction of the matter 
sources to the background can be ignored. This is the probe limit we will adopt in this 
paper. The background is taken to be a (3+ 1) dimensional Schwarzschild-AdS black hole 


ds? = — f (r)dt? + a + r?(dx* + dy’) (4) 
- f(r) Ši 
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with f(r) = r?(1 — h) and r, the horizon radius. The Hawking temperature for this 


black hole is T = Pira) = = which sets the temperature of the boundary field theory. 
In the probe limit, the matter fields p, and A, can be treated as perturbations on the 


Schwarzschild-AdS black hole background. 


3 Adding a constant magnetic field 


In Appendix A, we discussed the condensate of the vector field p, induced by the bulk 
electric field by turning on the scalar potential A; associated with the Maxwell field only. 
The normal phase corresponds to the black hole solution with a vanishing vector field p,.. 
As one lowers the temperature, the normal phase becomes unstable to developing non- 
trivial configuration of the vector field p,. It gives non-zero vacuum expectation value 
of the dual vector operator, which breaks the U(1) gauge symmetry and the rotational 
symmetry in x — y plane. The calculation of the optical conductivity reveals that there 
is a delta function at the origin for the real part of the conductivity, which means the 
condensed phase is indeed superconducting. For details, see Appendix A. 

We now turn on a magnetic field to study how the applied magnetic field influences on 
the system. A consistent ansatz is as follows ! 


podz” = [epz(r, x)” + O(e*)|dx + [ep,(r, zee” + O(e*)|dy, (5) 
A, dx” = [9(r) + O(€*)|dt + [Bx + O(e*)]dy, 
where pz(r, £), py(r, x) and ¢(r) are all real functions, p is a real constant, the constant 6 is 
the phase of p) and B > 0 is the constant magnetic field perpendicular to the x — y plane. 
We have defined the deviation parameter e€ from the critical point at which the condensate 
begins to appear. 
The zeroth order of (2) gives the equation of motion for @ 


g'e) +0) =0. (6) 


The asymptotic value of ¢ gives the chemical potential = A;(r — oo) of the dual field 
theory. The boundary condition at the horizon is given by requiring that A,A” is finite 
there. Thus we can obtain a unique solution 


Th 
a(r) = wa — 5). (7 
The equations of motion for py and py can be deduced from (3) at order O(e). We 
further separate the variables as p,(r,x) = Ya(r)X (x) and p,(r, x) = yy(r)Y (x). We find 
that to satisfy the equations of motion of the model with the given ansatz, 0 can only be 
‘One can also consider the dyonic black hole background with gauge field fixed [29]. p, is considered 


as a perturbation in such a background. In that case, the metric field f(r) will depend on electric charge 
and magnetic field B. 


chosen as 6, = 5 + 2nm or 0- = —5 + 2nz with n an arbitrary integer. The equations of 
motion for y:(r), yy(r), X(x) and Y (x) are divided into the following equations as 


PX + (qBx — p)p,Y =0, (8) 
pX + (qBx — p)p,Y =0, (9) 
fl ge? m2 : Y Y 
pat FP + Fy Pa — Pat raplF(aBe “Oy a +qBry — (gBx —p)’] = 0, (10) 
" a 1 Pp m? Py Y X Px X Yau 
L me it | (gBe —p)— 2 (1 B20, 11 
Py to Coe aa) peut aly (qBu YT a, (l+)q Yo, , (11) 


where the prime denotes the derivative with respect to r and the dot denotes the derivative 
with respect to x. Here and below the upper signs correspond to the 6, case and the lower 
to the 0_ case. In order to satisfy (8), one should impose 


Py = Chas X4 c(qBa — p)Y = 0, (12) 


where c is a real constant. This constraint is automatically satisfied by (9). We can see 
that only two of the four functions are independent. Substituting (12) into the remaining 
equations, we can find the following three equations 


f' eo? m2 E 
" / 
OA ot ee eee = (13) 
f P f r? f 
— X Fe(1+7)qBY + (Bz — p} X = EX, (14) 
ee penny Soe (15) 
Cc 


One can get the value of E for arbitrary constant c by solving the eigenvalue prob- 
lem (14) and (15) with the constraint given in (12). There may exist the possibility that 
one can not obtain non-trivial solutions for some special values of c. Here we consider a 
simple case with c? = 1, in which the equations of motion for X(x) and Y (a) can be solved 
exactly. Since c = 1/c, the c in the denominator in (15) is equivalent to the case in the 
numerator. Subtracting (14) from (15) and defining a new function as °? 


(a) = X(x) — Y (x), (16) 


one gets the equation 


p(s) + [E F c(1+7)gB — (qBz — p)"|p(x) = 0. (17) 


?To solve the above equations, one can also define y(x) = X(x) + Y (x). This case is equivalent to 
setting c + —c, which gives nothing new. 


: : EFc(1 B 
We further introduce a new variable € = \/|qB|(x — 4g) and a constant 7 = ——— 
then the above equation becomes 
2 


ev) + (n — EYE) = 0. (18) 
The regular and bounded solution of (18) is given by Hermite function H, as 
Wn lE) = Nne* /?H,,(€), (19) 


with the corresponding eigenvalue 7, = 2n +1. N, is a normalization constant and n is a 
non-negative integer. Thus we obtain the solution of equation (17) as 


Wn (2) = Nae 298-35)” H, (/|qBI (x — Bp) (20) 


with the corresponding eigenvalue 
E, = (2n+1)\q¢B\| + c(1 + y)qB. (21) 


Combining (12), (16) and (20), one can obtain the exact configurations for X(x) and 
Y (x), which read 


= eqB P 


cp2 
py =e"? (=~ a)" [X (0) e+ 208 + YBN x 


J e- Bye ae REE- a, 
— — je 2 qB T Aa ` 
0 qB ? qB 
and ; — 
Ya(2;p) = Xala: p) — Nae P-A? wy (aBa — 2), (23) 


qB 
where X (0) is a constant denoting the value of X,,(x;p) at the origin x = 0. The solutions 
of Yz and yy corresponding to the eigenvalue En, denoted by Yn and Yyn, can be obtained 
by solving the equation of motion (13) with Æ, given in (21). So far, we have recovered the 
Landau levels. As one can see in appendix A, we have a second order phase transition from 
the normal phase with p, = 0 to the condensed phase with p, # 0. Therefore we should 
encounter a marginally stable mode at the transition point. Theses solutions obtained 
from (13) just correspond to the marginally stable states. 
We can see from (13) that the effective mass of p, is 


2 2, Fn P _ 2 (n+1)lgB| +e +y)B Pg 
Z L io m G HUM L H, 
p2 f j2 J 
which is clearly shifted by the magnetic field B. Depending on concrete parameters and 
Landau level, the appearance of magnetic field can increase or decrease the effective mass, 


thus will hinder or enhance the transition from the normal phase to the condensed phase. 
In what follows, we consider the case with the lowest Landau level with n = 0, which reads 


Ey = —-|y4B], 


(24) 


No — 9B lg P_)2 (25) 


Xo (zip) =e 2 4B! = Yg (2; p). 


4 Phase diagram 


We are interested in how the applied magnetic field influences on the transition temperature 
from the normal phase to the condensed phase. As we can see from (24), the effective mass 
of the charged vector field p” in the lowest energy state, i.e., in the lowest Landau level 
n = 0 depends on the magnetic field B and the non-minimal coupling parameter y as 


ie =e = lyqB| E Pp 
eff r2 f 


It is clear that the increase of the magnetic field B decreases the effective mass and thus 
tends to raise the transition temperature. 


(26) 


4.1 Phase diagram at vanishing charge density 


It is well known that the increase of the electric field decreases the effective mass of charged 
scalar or vector fields, inducing the transition from the normal phase to the condensed 
phase. We first turn off the electric field, which corresponds to the case with vanish- 
ing charge density p = 0. We introduce a new coordinate z = r/r. The equation of 
motion (13) can be rewritten as ° 


n aa Oa eS o 
1—23"* z2(1— 23) 1672(1 — 2°) 


|~u(z) = 0, (27) 
with ¢ = |yqB|/T?. 

To solve such a second order equation by shooting method, we impose the regular 
condition at the horizon z = 1 as well as the source free condition at the boundary z = 0. 
More specifically, we set y,(1) = 1 in our numerical calculation due to the linearity of (27). 
For a given m?, only for certain values of ¢ = |yqB|/T? can the boundary conditions be 
satisfied. 

Figure 1 presents the three marginally stable curves of y,(z) for m? = 3/4. The three 
lowest-lying modes are in the sequence Co < ¢; < Gg. The red line corresponding to the 
minimal value of ¢ has no intersecting points with horizontal axis at non-vanishing z. Such 
a mode with Co ~ 72.84 is considered as a mode of node n = 0. Furthermore, the green 
line to ¢, œ 339.24 and blue line to ¢2 ~ 780.10 are regarded as modes with nodes n = 1 
and n = 2, respectively. Since the radial oscillations in z-direction of y,(z) will cost more 
energy, the later two curves are therefore thought to be unstable. Thus the lowest value Co 
just gives the critical magnetic field above which the normal state is unstable to developing 
a vector hair. Figure 2 shows the critical magnetic field in terms of Ço for various squared 
mass of the vector field. It can be seen clearly that Co increases as we increase the squared 
mass. 

3We apologize to the readers for here using a same notation to denote the derivative with respect 


to different variables for brevity. But the meaning of the derivative in the text is clear and will not be 
confused. 


Z 


Figure 1: The marginally stable curves of the charged vector field corresponding to various 
critical ¢ = |yqB|/T?. The three curves from top to down correspond to Go ~ 72.84 (red), 
Cı œ 339.24 (green) and 2 ~ 780.10 (blue), respectively. We choose m? = 3/4. 
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Figure 2: The critical magnetic field Çọ versus m^ of the vector field. The points are 
obtained by the shooting method to solve the equation (27). 


It should be pointed out that if we turn off the magnetic field, the normal state will 
not become unstable to developing hairs with vanishing charge density. The interesting 
result here is that only the magnetic field itself can trigger the phase transition. This 
result has an analogy to the QCD vacuum instability induced by a strong magnetic field 
to spontaneously developing the p-meson condensate. It is clear that the last term in (1) 
describing a non-minimal coupling of the vector field p“ to the gauge field A, plays a crucial 
role in the instability. Note that similar coupling can be found in many formalisms used 
to describe the coupling of magnetic moment to the background magnetic field for charged 
particles of spin 1, i.e., vector particles [19, 32]. 


10 


4.2 Phase diagram at finite charge density 


Let us now consider the system with fixed charge density p. Here we do not limit ourselves 
to the lowest energy state. The equation of motion (13) can be written in terms of p and 
Zz as 
327 m? (En/ PA gr? 
" / n 
(2) = = | ay ee le) =; 


(28) 


where A= 4 and En, is given in (21). The lowest energy state corresponds to En = Spe 


—|yqB|. For numerical convenience and to match the behavior at the boundary, we further 
define 


palz) = (—)*-F (2). (29) 
Then we can obtain 
PM) = 2 2) Ig BS - Sr 
Ien (30) 
E O EDA o 


4.5p 
4.0 F 
3.56 normal phase 
3.0} 
2.5 F 
2.0} 
1.5} 
1.0} 


T/T. 
condensed phase 


lyqB/p| 


Figure 3: The transition temperature from the normal phase to the condensed phase as a 
function of magnetic field. It corresponds to the case with En = —|yqB]. Te is the critical 
temperature in the case without magnetic field. The magnetic field raises the transition 
temperature. We choose m? = 3/4. 


This equation depends on two physical parameters ie = ink and a To solve such 
second order equation, we impose the regular condition at the horizon z = 1 as well as 
the source free condition F(0) = 0 at the boundary z = 0. It has a non-trivial solution 
only when there is a relation between such two parameters, which just gives the transition 


temperature as a function of the magnetic field. 
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The (T, B) phase diagram with the lowest Landau level is drawn in figure 3. To deter- 
mine which side of the phase transition line is the condensed phase, we can consider the 
equation (26). It suggests that the magnetic field decreases the effective mass. So if we in- 
crease the magnetic field at a fixed temperature, the normal state will become unstable for 
sufficiently large magnetic field. Figure 3 looks very similar to figure 9 in ref. [22] where the 
chiral transition temperature rises with magnetic field, indicating chiral magnetic catalysis. 
Furthermore, to compare with the lowest Landau level case, we present an example with 
positive En = |yqB| in (28) in figure 4. One can see clearly in this case that the transition 
temperature lowers with the increase of the applied magnetic field. It is the well known 
property of the ordinary superconductor, which has been first discussed in a holographic 
setup in ref. [33]. 


1.0F 


0.97 


Der normal ph 
T/T. ormal phase 


0.7} 
0.6} ] 
condensed phase 


0.5 El 1 1 in L 1 į 1 1 1 Í 1 L 1 1 1 L ï 1 1 a 
0 2 4 6 8 10 


lyqB/p| 


Figure 4: The transition temperature from the normal phase to the condensed phase versus 
magnetic field. It corresponds to the case with En = |yqB|. Te is the critical temperature 
in the case without magnetic field. The magnetic field leads to lowering the transition 
temperature. We choose m? = 3/4. 


5 Vortex lattice solution 


Let us now construct the vortex lattice solution. It is enough to consider the n = 0 solution 
only, i.e., 
1 2 
polz; p) = Noe 2/9710 aa” (31) 


Due to the fact that the eigenvalue E,, is independent of p, a linear superposition of the 
solutions e”9 pen(r)Xn(2; p) and ec”! pyn(7)Y;,(x; p) with different p is also a solution of the 
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model at O(e). We introduce two functions 


+00 

Pelr, £, Y) = Pool p> ce XG (x; pe), 
L=—0o 
+00 

Annum ee aalr DD ce Oy (oP); (32) 

l=—0o 

4292 92 27 /\qBlé 

ce =e , Pe= a | , 

ay 


which satisfy the full equations of motion. a, and as are arbitrary constants. Following 
ref. [33], we can obtain the vortex lattice solution 


+00 


*py(r, 2, Yy) = Pxo(r) » cee"? Wo (x; pe). (33) 


L=—00 


—i0 


Ap(r, x,y) = palr, £, y) — ce 


Ap(r,x,y) has a pseudo-periodicity 


Ap(r, 2, y) = Ap(r, £, y + ), 


Toe 


2 idm ae (34) 
Apr, 2 + ——, y + A) = est VBA Apr, 2, y), 
ar y qB] aıy/lqB| 
as well as a zero at 
Xma = (m s)bi + (n+ Ib», (35) 
with two vectors b; = Jan” and b aia" Faa m and n are two integers. 


Since the expectation value of the operator J” dual to Pp is given by the coefficient 
at boundary r > oo, the quantity J+ = (Je + if”) indeed exhibits the vortex structure 
with the cores of vortices located at Xm,n. In particular, the triangular lattice with three 
adjoining vortices forming an equilateral triangle can be obtained by choosing the following 


parameters 
a, = BLN io = all (36) 
V3 V3 
It should be stressed that it is the special combinations J+ which exhibit the vortex lattice 
structure. In particular, for q > 0, it is J. that corresponds to the lowest Landau level, 
while for q < 0, it is J}. The form of operator presenting vortex lattice structure is the 
same as the one in field theory study without gravity [16]. This also provides the evidence 
for the correctness of choosing c? = 1 in section 3. 
Figure 5 shows the configuration of the norm of condensate J_ in the x — y plane for 
the triangular lattice. Obviously, to obtain the true ground state, we should calculate the 
free energy of the solutions with different lattice structures from the action to find which 
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configuration minimizes the free energy. It turns out that the linear analysis presented 
here is not sufficient to determine the most stable solution. We should include higher order 
contributions just as done in refs. [24, 33]. The calculation is much more complicated and 
is not very relevant to our purpose of this paper. We leave it for our further study. 


Figure 5: The vortex lattice structure for the triangular lattice in x—y plane. The contour 
plot is also drawn in the bottom. In particular, the condensate vanishes in the core of each 
vortex. 


6 Conclusion and discussion 


In this paper we studied a holographic model with a complex vector field charged under 
a U(1) gauge field in a (3 + 1) dimensional AdS black hole background, aiming to shed 
some light on the real strongly coupled systems which are of gravity duals. In this model, 
there is a non-minimal coupling of the vector field to the U(1) gauge field, which describes 
the interaction between the magnetic dipole moment of the vector field to the background 
magnetic field. This model includes the minimal ingredients to build a holographic p-wave 
superconductor model. We found a critical temperature at which the system undergoes a 
second order phase transition. The critical exponent of this transition is one half which 
coincides with the case in the Landau-Ginzburg theory. In the condensed phase, a vector 
operator acquires a vacuum expectation value breaking the Abelian symmetry as well as 
rotational symmetry spontaneously. Our calculation showed that this condensed phase 
exhibits an infinite DC conductivity and a gap in the frequency-dependent conductivity, 
which is quite similar to properties of the ordinary superconductivity. 
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We paid more attention on the response of this system to an applied magnetic field. We 
obtained the Landau level, from which we can find the contribution to the effective mass of 
the vector field by the magnetic field (see (24)). Due to the non-minimal coupling given in 
the last term of (1), the applied magnetic field can reduce the effective mass of the vector 
field, thus inducing the instability of the black hole background even when the chemical 
potential /charged density is absent. That is, for the case with vanishing chemical potential 
or charged density, the black hole background becomes no longer stable when the magnetic 
field is beyond a certain critical value. For the case with non-vanishing chemical potential, 
the phase boundary is determined by a relation between the transition temperature and 
the magnetic field, which is presented in figure 3. The transition temperature increases 
with the applied magnetic field. The response of this system to the magnetic field is quite 
different from the behavior of ordinary superconductor where the magnetic field makes the 
transition more difficult as drawn in figure 4. But our result is quite similar to the case of 
QCD vacuum instability induced by strong magnetic field to spontaneously developing the 
p-meson condensate [16, 17]. Although so, it was shown that in our model, the condensate 
of the vector operator forms a vortex lattice structure in the spatial directions perpendicular 
to the magnetic field. Of course, the non-minimal coupling term in the action plays a crucial 
role in both cases. Therefore in some sense, our model is a holographic setup of the study 
of p-meson condensate in refs.[16, 17]. 

In ordinary superconductors an external magnetic field suppresses superconductivity 
via diamagnetic and Pauli pair breaking effects. However, it has also been proposed that 
the magnetic field induced superconductivity can also be realized in type-II supercon- 
ductors [34, 35], in which the Abrikosov flux lattice may enter a quantum limit of the low 
Landau level dominance with a spin-triplet pairing. And possible experimental evidence for 
the strong magnetic induced superconductivity can be found, for example, in refs. [30, 31]. 
It was also shown in Gross-Neveu type model that applied magnetic field might induce 
superconductivity in the planar system with 4-fermion interaction [36]. 

We mention here that similar studies can also be generalized other gravitational back- 
grounds, such as the AdS soliton background which has been adopted to mimic supercon- 
ductor/insulator phase transition [37]. The study of magnetic field effect in the supercon- 
ductor/insulator case can be found in ref. [38]. In a forthcoming paper [39], we generalize 
the present study to the case with the AdS soliton as the background. It is found that the 
magnetic field can induce the instabilities of the AdS soliton background. Comparing our 
model with the SU(2) model with a constant non-Abelian magnetic field [23, 25], we find 
that our complex vector field model in some sense is a generalization of the SU(2) model 
to the case with a general mass squared m? and magnetic moment characterized by y. In 
the setup of the present paper, the SU(2) model corresponds to our model with m? = 0 
and y = 1. 

In this paper, we restricted ourselves to the probe approximation, neglecting the effect 
of matter fields on the background geometry. This can indeed reveal some significant prop- 
erties of the model, but something might be lost in this approximation, see ref. [40] as an 
example. It is therefore helpful to understand full properties of the model by considering 
the back reaction of matter fields on the background geometry. In a recent paper [41], going 
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beyond the probe approximation, we found a rich phase structure in this model without 
magnetic field. Depending on mass square m? and charge q of the vector field, not only 
second order but also first order and zeroth order phase transition can appear. Interest- 
ingly, there also exists a so-called “retrograde condensation” in which the hairy solution 
exists only for temperatures above a critical value and is thermodynamically subdominant. 
Particularly, the zeroth order transition and retrograde condensation can be observed for 
the case with small m?. Indeed this model has much more phase behaviors than the SU(2) 
model, thus it can be used to mimic much richer phenomena in dual strongly coupled 
systems. 

In our model, it is the magnetic field itself that can induce the condensate. A natural 
question arises how about the Meissner effect known as that superconductors expel weak 
external magnetic field. This effect is due to the large superconducting currents induced 
in the superconductor by the external magnetic field, which generates a back-reacting 
magnetic field screening the external magnetic field. In the holographic SU(2) p-wave 
model, it has been pointed out [25] that the vortex currents flow in the opposite direction 
to the one in conventional superconductors, thus can enhance the applied magnetic field in 
the regions between the vortices. It seems that the model presented here is similar to the 
SU(2) p-wave model. It is interesting to ask whether a similar phenomenon also appears in 
our model. To answer this question, we need to consider the contribution of higher order 
terms in (5), which is also required to find the true vortex flux lattice. We leave it for 
further study. 

In this paper we only calculated the conductivity in the case without magnetic field in 
one spatial direction. Although it is enough to see the superconductivity feature of the 
condensed phase, in order to study the model in a more realistic manner, it is desirable to 
calculate the conductivity in another direction and further to study the transport properties 
of the lattice state. As a phenomenological approach, this toy model would have potential 
application to mimic strongly coupled systems with a vector like order parameter. It should 
be applicable in a wide variety of condensed matter systems, heavy ion physics and beyond. 
We hope to report further progresses in future. 
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A Condensate of vector field and superconducting phase 
transition 


As a toy model with a charged vector field in the bulk dual to a vector operator in the field 
theory, there exists the possibility that the condensate of this vector field can serve as an 
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order parameter to mimic a holographic p-wave superconductor phase transition, like the 
s-wave case |4, 5]. More precisely, we hope that this system has stable black hole solutions 
with vector hair at low temperatures, but without vector hair at high temperatures. If 
this is true, in the condensed phase, the condensate of the dual vector operator will break 
not only the U(1) symmetry but also the rotational symmetry since the condensate of 
vector field picks out one special direction. This situation is very similar to the one in the 
holographic p-wave superconductor model with SU(2) gauge field in ref. [12]. In this sense, 
our model can also be regarded as a holographic p-wave model. 


A.1 Condensate of vector field 
We adopt the following ansatz 


pdz” = pz(r)dx+ py(r)dy, A dz” = d(r)dt. (37) 


One can use the U(1) gauge symmetry to set py to be real. Then one finds that the r 
component of (2) implies that the phase of p, must be constant. Without loss of generality, 
we take p, to be real. Then, the independent equations of motion in terms of the above 
ansatz are deduced as follows 


/ 2 2 
o! TETEE Cpe MPs A, 
! 242 2 
g 
p Ty- 7 (02 + #4)0 = 0, 


where the prime denotes the derivative with respect to r. 

In order to find the solutions for all the three functions F = {p,, p,,} one must impose 
suitable boundary conditions at the AdS boundary r — oo and at the horizon r = rp. In 
addition to f(ra) = 0, one must require ¢(r;,) = 0 in order for g” A,A, to be finite at the 
horizon. 

In order to match the asymptotical AdS boundary, the general falloff of the matter 
fields near the boundary r — oo should behave as 


aya! — Pie, Per = Pra y e 
p =ü sis hm a a M= ta (39) 


where Ay = oat ts 4 We impose p,_ = 0 and py_ = 0, since we want the U(1) 


symmetry to be broken spontaneously. According to the AdS / CFT dictionary, up to a 
normalization, the coefficients 4, p, Pr} and py, are interpreted as chemical potential, 
charge density and the x and y components of the vacuum expectation value of the vector 
operator J! in the dual field theory, respectively. 


4The m? has a lower bound as m? = —1/4 with Ay = A_ = 1/2, in which case there is a logarithmic 
term in the asymptotical expansion. Such a term should be considered as the source set to be vanishing [42]. 
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There is a useful scaling symmetry in the equations of motion 


r= Ar, (t,2,y) >A *(t,2,y), ($, Pes Py) > MO, Px, Py); (40) 


where A is an arbitrary positive constant. Under this symmetry, the revelent quantities 
transform as 


TOM, pr, pp, (Pap Pyy) AAMT Pr4 Py): (41) 
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Figure 6: The condensate as a function of temperature. We choose q = 1 and m? = 3/4. 
The condensate begins to appear at Te ~ 0.102,/p. 


We assume the condensate to pick out the x direction as special, so we can consistently 
set py = 0. The condensate as a function of temperature is presented in figure 6. It is clear 
that as we lower the temperature, the normal phase with vanishing charged vector becomes 
unstable to developing vector hair which breaks the U(1) symmetry as well as rotational 
symmetry in the dual field theory. Fitting the curve near the critical temperature Te, we 
find that for small condensate there is a critical behavior with critical exponent 1/2, which 
precisely meets the result given by mean field theory and it is typically a second order 
phase transition. In the case with m? = 3/4 we obtain 


(i) = 339T5P (1 — T/T.) P, as T >T.. (42) 


Thus we have obtained two black hole solutions in the system. When T > Te, we 
have the black hole solution without the vector field, while we have the black hole solution 
with non-trivial vector field as T < Te. This behavior is in complete agreement with 
the holographic superconducting phase transition in the literature. Therefore we expect 
that the black hole solution with non-trivial vector field can describe a superconducting 
phase. To prove this, it is helpful to calculate the optical conductivity. Before dong this, 
we should first ensure that the black hole solution with non-trivial vector field is more 
thermodynamically stable than the one without the vector field as T < Te. 
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A.2 Free energy 


In order to determine which solution is thermodynamically favored, we should calculate the 
free energy of the system for both black hole solutions. We will work in canonical ensemble 
in this paper, where the charge density is fixed. In gauge/gravity duality Helmholtz free 
energy F of the boundary thermal state is identified with temperature times the on-shell 
bulk action with Euclidean signature. Since we work in the probe approximation, we can 
ignore the gravity part. Given that the system is stationary, the Euclidean action is related 
to the Minkowski case by a total minus. Employing the equations of motion (2) and (3), 
we have 


—2* S Euclidean = [ar —g( 


YV 1 V s. V 
Fgh = P =m pie Figer] 


ii 
4 2 


+ I dry —hn ALF + Sa (43) 


1 1 
= | PaE AN PO + perv Thin (Auk — pi pt”) + Sct, 


where h is the determinant of the induced metric h,, on the boundary r —> oo and n“ is 
the outward pointing unit normal vector to the boundary. Se denotes the surface counter 
term for removing divergence. 

Substituting the asymptotically expansion (39) into (43) and introducing a counter 
term Sa = —A_ f dxz?V—hh,vp""p”, ° we find the free energy F as 


27°F 1 i 1 mA 
= 5HP 7T (2A, a 1)(pPz-Pr+ + Py-Py+) = dry 9, AV LF , (44) 
V 2 : 2 


h 


with V = f drdy. Regarding p,— and p,— as sources, Helmholtz free energy tells us that 
sub-leading terms p,, and p,, are the expectation values of the dual operator in the field 
theory side. The difference of Helmholtz free energy between the condensed phase and 
the normal phase as a function of temperature is presented in figure 7. It is clear that 
below the critical temperature T}, the state with non-vanishing vector “hair” is indeed 
thermodynamically favored, compared to the normal phase. The phase transition is second 
order, which can be seen, for example, from the derivative of the free energy with respect 
to the temperature. 


A.3 Conductivity 


We now calculate the conductivity in the dual field theory side as a function of frequency 
w. We need to turn on fluctuations of the matter contents in the bulk. We assume 
perturbations have a time dependence of the form e~*? with zero spatial momentum. It 


>We are not sure whether the counter term works or not in a general case. But we find this counter 
term works well for the ansatz (37) with the asymptotically expansion (39). 
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Figure 7: The difference of Helmholtz free energy between in the condensed phase and in 


the normal phase as a function of temperature. We choose q = 1 and m? = 3/4 for which 
the critical temperature Te ~ 0.102,/p. 


turns out that one can consistently turn on the perturbations of A, only. We can obtain 
the equation of motion for A, by linearizing the equations of motion (2), which reads 

P w2 oy a 

Ad / T 

aT (z7 rf 

Since the conductivity is related to the retarded two-point function of the U(1) current, 

we impose the ingoing boundary condition near the horizon. The gauge field A, near the 

boundary r — oo falls off as 


JAn (45) 


AW 
A, = A + —+---. (46) 
r 


According to the AdS/CFT dictionary, one can obtain the conductivity as 


AW 
Oyy(W) = iw A) ` (47) 


The AC conductivity as a function of frequency is presented in figure 8. We can see 
clearly that the optical conductivity along the y direction in this model behaves qual- 
itatively similar to the case in the p-wave model with SU(2) gauge symmetry [12]. In 
particular, from the Kramers-Kronig relation, one can conclude that the real part of the 
conductivity has a Dirac delta function at w = 0 since the imaginary part has a pole, i.e., 
Im|[oyy(w)] ~ +. Furthermore, it is clear that the optical conductivity develops a gap at 
some special frequency w, known as gap frequency. As suggested in ref. [42], it can be iden- 
tified with the one at the minimum of the imaginary part of the AC conductivity. Re[a,,| is 
very small in the infrared and rises quickly at wg. There also exists a small “bump” slightly 
above w, which is reminiscent of the behavior due to fermionic pairing [12]. For our chosen 
parameter, we have wg ~ 8T.. Compared to the corresponding BCS value wg ~ 3.57%, the 
result shown here is consistent with the fact that our holographic model describes a system 
at strong coupling. 
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Figure 8: The optical conductivity as a function of frequence. The solid lines in the left 
plot are the real part of the conductivity, while the dashed lines in the right plot are 
the imaginary part of the conductivity. We choose q = 1 and m? = 3/4. The horizontal 
lines correspond to the temperature above T}. Other curves from left to right correspond to 
T/T, ~ 0.830 (purple), T/T, ~ 0.519 (green), T/T, ~ 0.388 (blue), and T/T, ~ 0.290 (red), 
respectively. There is a delta function at the origin for the real part of the conductivity in 
the condensed phase. 
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